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Abstract 

Elastic electromagnetic form factors of nucleons are investigated both for the time-like 
and the space-like momentums by using the unsubtracted dispersion relation with QCD con- 
straints. It is shown that the calculated form factors reproduce the experimental data reason- 
ably well; they agree with recent experimental data for the neutron magnetic form factors for 
the space-like data obtained by the CLAS collaboration and are compatible with the ratio of 
the electric and magnetic form factors for the time-like momentum obtained by the BABAR 
collaboration. X 

PACS No. 13.40. Gp, 12.38.-tni, 11.55.Fv 

1 Introduction 

Recently, there are remarkable developments in experiments for the nucleon electromagnetic form 
factors: 

1) For the space- like momentum, the ratio of the electric and magnetic form factors of pro- 
ton, and G'^j respectively, was shown to be a decreasing function of the squared momentum 
transfer and the experimental results imply that the proton electric form factor vanishes for 

« 7(GeV/c)2 

2) For the neutron magnetic form factor, G^^, very accurate experimental data were obtained 
and it approximately satisfies G^,^(Q^)//z„ ~ Gd{Q^) = (1-I-Q^/0.71)~^, with Q being represented 
in terms of GeV/c, for fairly wide range of squared momentum transfer = 1.4 — 4.8 (GeV/c)^ 
[g, [7] (CLAS collaboration). 
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3) For the time-like momentum the ratio jG^/G^^I was obtained [5], [5] (BABAR cohabora- 
tion), while previously the data of form factors had been analyzed under the assumption = 
or G^ = G^f . 

Asymptotically, the experimental data of nucleon magnetic factors decrease more rapidly than 
the dipole formula for large and the decrease has been understood as a realization of pertur- 
bative QCD [TO] , the behavior of which can be formulated in terms of the dispersion theory with 
appropriate conditions on the absorptive parts; we assumed unsbtracted dispersion relations for 
the charge and magnetic moment form factors. To realize the asymptotic form of QCD we imposed 
the superconvergence conditions. 

As the data for the time-like momentum have become accurate, it is necessary to investigate 
the form factors for the space-like and time-like momentums systematically. For this purpose the 
dispersion theory is effective. 

The dispersion theoretical calculations performed so far, the value of G\i turned out to be 
larger than the above mentioned new experimental data for = 1.4 — 4.8 (GeV/c)^ (see Ref. [7]). 
It is of vital importance to investigate if it is possible to realize the experimental data simply by 
the adjustment of parameters or by the refinement of absorptive parts in the dispersion relation. 

It is the purpose of this paper to analyze experimental data of nucleon form factors by the 
dispersion theory, with the QCD constraints imposed, taking account of the above mentioned new 
experimental results. 

Organization of the paper is given as follows: In Sec. 2 we explain the superconvergent dispersion 
relation and give conditions which are used in this paper. We summarize the absorptive parts, 
which are broken up into three parts: Low, intermediate and asymptotic momentum regions. 
For each momentum region the imaginary parts are given. The asymptotic part is expressed as an 
expansion in terms of the analytically regularized running coupling constant in the renormalization 
group for QCD. In Sec. 3 we remark on the numerical analysis. In Sec. 4 numerical results are 
summarized. The final section is devoted to general discussions. 

2 Dispersion Relation for the Electromagnetic Form Factors 

We assume the unsubtracted dispersion relations for the charge and magnetic moment form factors, 
FI and , respectively, with / denoting the isospin state / = 0, 1. That is. 



where the threshold is to — 4/i^. Here /i is the pion mass being taken as the average of the neutral 
and charged pion masses. We impose conditions on ImF/ to realize the QCD conditions. 

2.1 Superconvergence Condition and QCD 

Experimental data imply that the magnetic form factors of nucleon decrease more rapidly than 
the dipole formula for large squared momentum transfer. The decrease agrees with the prediction 
of perturbative QCD, where magnetic form factors of nucleon decrease for — > oo as 



where ag is the running coupling constant of QCD and /?o — 11 — 2n//3 with rif being the number 
of flavor. A is the QCD scale parameter having the dimension of momentum. 




(1) 




(2) 
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To realize the QCD predictions we impose the fohowing conditions on the charge and magnetic 
moment form factors: 

Fi(Q2) ^ const/[Q2(lnQ7A2)^], 

F2iQ^) -> const/ [Q4(lngVA2)^], (3) 

for — 7^ CXI with 7 > 2. 

We briefly summarize the asymptotic theorems which are used to incorporate the constraints 
of QCD [TU], where the proof is given in Ref. [TB]. Let F{t) satisfy the dispersion relation (IT|), and 
Imi^ is given as 

= R^V^^ + ^ ( [ln(tvl^)]7+0 
for t ^ 00 with 7 > 1. Then F{t) becomes 

1 r°° c 

c 



7r7 ln(|t|/A2)]7 
for t — > ±00. Generally, when F{t') satisfies 

for — )■ 00 and the superconvergence conditions 

I dt't'^lmF{t') = 0, fc = 0, 1, • • • , n, (7) 

J to 

F{t) given by ([T]) approaches for t ±00 to the following formula: 



^^Jto t'-t i"+i^7[ln(|t|/A2)7' 



which can be proved by using ^ and ([7]) together with the identity 

= --<! 1 + 

t' 

Indeed, by using ([7]) we have 



' dt''^ ^ ' r (9) 

t' - i t"+i ./to t' -t ' ^ ^ 



which leads to © as t'"+ilmF(t') satisfies dH). 

To obtain the asymptotic formulas ([3]), therefore, we impose the superconvergence conditions 
on the imaginary part of form factors, ImF^{t) [i = 1, 2; / denotes isospin) in the unsubtracted 
dispersion relation ([T]): 

-, /.oo A r°° 

- dt'lmFl{t') = -/ dt't'lTnFl{t')=0, (10) 

J to ^ J to 

-/ dt'lmFiit') ^ - dt't'lmFiit') ^ - dt't'^ ImF^ {t') ^ 0, 

J to J to ^ J to 
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where ImF/(i') satisfies the asymptotic conditions for t' ^ oo 

t'nuiF^it') const/[ln(f7A2)]''+i {i = 1, 2). (11) 
In addition to the conditions ([TU]) and ((TT|) we impose the normalization conditions at t = 0: 

- = - / dt'lu,Fi{t')/t', (12) 

^ to 

1 

g' = - dt'lmFi{t')/t\ (13) 

^Jto 

where is the anomalous magnetic moments of nucleons with the isospin /. 
2.2 Imaginary part of the form factors 

Let us discuss the imaginary parts of nucleon form factors, which are broken up into three parts: 
The low momentum, the intermediate, and the asymptotic regions. 



2.2.1 Low momentum region 

The imaginary parts of the charge and magnetic moment form factors, ImF^^, are given in terms 
of two pion contribution as follows: 



lm[Fr{t)/e] = 



m (t- V) /^i-4^2\ 1/2 
"2 4m2 - t 



X Re 



M 



t m (_) 



Am? V2 



lm[2mF^{t)/e] 



2 (4m2 -t) \ t 



(14) 



X Re 



IV2 



where /j_ ^^{t) are helicty amplitudes for tttt o NN, M(t) is the pion form factor and /i is the 
pion mass. The superscript V denotes the iso-vector part. For the helicity amplitudes we use the 
numerical values given by Hohler and Schopper |16) and parameterize M{t) according to them. 



M{t) = tp{l + {Tp/mpd)}[tp - t - tmlTp{qt/qpfVt]-\ 
where rUp and Tp are the p meson mass and width respectively and 

tn 



rrip, Qp 



tp - M^, 



d = 



TTt„ 



In- 



2qp 



2fi 



2Tiq, 



(15) 

(16) 
(17) 



The imaginary parts thus obtained are denoted as ImF/^ (j = 1, 2) hereafter. It must be remarked 
that the p meson contribution is included in the helicity amplitudes of Ref. [TB] . The uncorrelated 
kaon pair is neglected here as the effect was estimated to be small [TT] . 
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2.2.2 Intermediate region 

The intermediate states Afi^ < t < are approximated by the addition of the Breit-Wigner terms, 
with the imaginary part parameterized as follow: 

where 

rM|(M| + lit -tor i 



Here Mji and F are the mass and width of resonance, respectively, the threshold to is to = and 
tres is treated as an adjustable parameter, g is introduced to cut-off the Breit-Wigner formula. 
We write the intermediate part as the summation of resonances 

hnFr^' ^J2^'fnR' (20) 

n 

where / is the isospin and n is the labeling of resonances (see Table I). Here the sufSx i denotes 
i = 1, 2, corresponding to the charge and magnetic moment form factors and (N = n or 
p). The same formulas for /^^ are used for z = 1 and i — 2. 

2.2.3 Asymptotic region 

We express the form factors as power series in the running coupling constant of QCD, as- To 
calculate the absorptive part, it is necessary to perform analytic continuation to the time- like 
momentum. Here we give only the necessary procedure for the analytic continuation of the running 
coupling constant to the time-like momentum by using the analytic regularization |12j |13| , as the 
formulation is given in Ref. [18 . 

Let as{Q^) be the running coupling constant in the renormalization group calculated by the 
perturbative QCD as the function of the squared momentum for the space-like momentum. We 
use the three loop approximation for as{Q'^), which is expressed in the Fade form. 

MQ-') ^ I [>..«3VA^) + «■ Infln(QVA^)} + ".^tl^ + 5^ + ^ ^ ^ 

A is the QCD scale parameter, and are expressed in terms of the /3 function of QCD, 

where 



al = 2/3l//3o^ 02 = 4^, «3 = ^ ( 1 - ^ ) , (22) 



/5o = ll-^, /3.=51-l|^,/3.. 2357-^n, + f (23) 

with Uf being the number of flavor. We perform the analytic continuation of the squared momen- 
tum to the time- like region, s, by the replacement in (|2ip 

Q2 ^ e-'^'s. (24) 

Then as{e~^'^s) becomes complex and is expressed as follows: 

II -\- 111 

asie'^^'s) = l/{u-tv) = —j^, (25) 
D = u^+v^, (26) 
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where u and v are given as 



ln(s/A2) + ^ In{ln2(s/A2) + tt^} 



02 



In2(s/A2) +7r2 

a3ln(s/A2) 
"ln2(s/A2) + 7r2 



In(s/A2)ln{ln2(s/A2) +7r6l} 



(27) 



02 



ln^(s/A2) +7r2 

TTOs 

ln{ln2(s/A2) + 7r2}' 



2 In{ln2(s/A2) + TT^} - 01n(s/A2) 



(28) 



with 

61 = tan-i{7r/ln(s/A2)}. (29) 

The running couphng constant is given by the dispersion integral both for the space-hke and the 
time-hke momentum 

aa{t)= dt'^^ (30) 

with 

cr(i') = Imas(e~"s) = ^ttv/PoD. (31) 

aii{t) represented by (j30p is called analytically regularized running coupling constant as it has no 
singular point for t = — Q2 < Q. The regularization eliminates the ghost pole of as(Q^), given by 
(HD), appearing at 

Q2^g*2^^2g«*^ (32) 

where u* — 0.7659596 • • • for the number of flavor Uf = 3. Calculating pOL we find that aii{t) is 
approximately given by the simple formula with the ghost pole subtracted 

aRm^as{Q^)-A*l{Q^-Q*^), (33) 

where the residue A* is 



...4,AV/{*(l + ::^-e,i^ + ^)}^ (34, 



We use ((33l) as the regularized coupling constant; for the time-like momentum we replace 
e~*^s in ([55)) as was mentioned before. 

The QCD parts, F^'^^'^ {i — 1,2; I = 0,1) for the squared time-like momentum, are written 
as follows: 

j,QCA/(,)^^QCA/(,);,^(,)^ (35) 

where Pf^'^^'^'s are given as expansion in terms of the running coupling constant 

J>2 
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for the time-like squared momentum s. We multiply by the function h{s) in ([55]) to assure the 
convergence of the superconvergence conditions (llOp and (|lip . The following formula is assumed 
for hi{s): 

, X 3/2 / , s i+1 



which may be interpreted as the form factor for 7 — >■ with tg being the threshold of the quark 
antiquark pair. The parameters tq, ti and ^2 are taken as adjustable parameters and will be 
determined by the analysis of experimental data. 

For the time-like momentum, we perform the analytic continuation of the regularized effective 
coupling constant afi{Q^) to aii{s) through the equation 

anis) = afl(Q2e— ) = Re[aR{s)] + tlm[aR{s)]. (38) 

We express the QCD part as the power series expansion in ck_R(s) 

Fr"^'is) = Y.€^^'Ms)y- (39) 

The summation in begins in the second order in the effective coupling constant so as to realize 
the logarithmic decrease of the nucleon form factors. 
Imaginary part of (P^]) is obtained to be 

ImF^ ' — 2c^^^'^Reafl,Ima/j 



-c. 



QCD, I 
'i,3 



+cf^i°'^[4{ReaRflmaR - 4Reafl(Im aii)^] 

+ ••'•, (40) 

and 

Imi^'3^^'^(s) = lmF^'^^^\s)h,{s). (41) 

We write the low energy part, intermediate resonance part and asymptotic QCD parts of form 
factors as F^, pf^'^ and p^^^^^ ^ respectively, which are given by the dispersion integral with 
the imaginary parts p^ . ([20)1 and (HI]). The form factors F/ are defined by adding them up. We 
impose the conditions (ITUl) and (ITT|) on ImF/ so that the QCD conditions are satisfied. 



3 Numerical Analysis 

We analyzed the experimental data of nucleon electromagnetic form factors G\j/ ^.pGo, G'^/Gd^ 
G%[/ i^nGo G'% and the ratio iipG^/G^j for the space-like momentum transfer, and jCj and jG"! in 
Refs. |19) - |40j for the time-like momentum transfer and the above mentioned recent experimental 
data for the space-ike and |G^//ipG^^| for the time-like momentum transfer. The parameters 
appearing in the formulas are determined so as to minimize x^. 

As was mentioned in the introduction we analyze by taking account of the recent experimental 
data: (a) Gli for ^ 1 - 4.8 (GeV/c)^ (CLAS collaboration) and (b) |AipG|l/|G^,f I (BABAR 
collaboration) . 

In order to see how the situation changes by taking account of these new experiments in addi- 
tion to the other data, we perform analysis for the following two cases in the analysis: 
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Case I: Both of the experimental data, (a) j/ipG^/G^j^l for the time-like momentum and (b) 
new data for G'^j for the space-like momentum, are added. 

Case II: Only the data (a) j/ipG^/G^^I for the time-like momentum are added. 

Let us remark on the experiments for the time-like momentum [S], |32| . |33| . where the form 
factors \G'P\ and |G"| are determined by using the formula for the cross section cjo for the processes 
e-fe— ^7V-|-7VoriV + 7V— J-e-l-e, which is given as 



2 



3s 



1 + — |G(s)p. (42) 



Here a is the fine structure constant and v is the nucleon velocity. |GjJ^| are estimated from |G| 
under the assumption Gm = G^ or Gb = 0. do is expressed in terms of Gj^ and as follows: 

.o = ^(|G-p + ^|Gf|^). (43) 

Equating and (gS]), we have 

, _ |Gf,p + 2m^|Gg|V. 
- l + 2mVs • ^^^> 

Substituting our calculated result of form factors to the right hand side of (H^ . we obtain the 
theoretical value for |G|, which is compared with the experimental data for the magnetic form 
factor obtained under the assumption Gm — Ge- 

The parameters appearing in our analysis are the following: Residues at resonances, coefficients 
appearing in the expansion by the QCD effective coupling constants, cut-offs for the intermedi- 
ate region Ai. In addition to them we have parameters in the Breit-Wigner formula and the 
convergence factor h of QCD contribution, tg, t^es, ii, ^2, ^3- 

We have taken the masses and the widths of resonances as adjustable parameters. As the su- 
perconvergence constraints impose very stringent conditions on the form factors, it was necessary 
to take the masses and widths as parameters. 



4 Numerical Results 

We give in Tables 1, 2 and 3 the results for the parameters for the cases I and II obtained by the 
analysis; in Table 1 the masses and widths of resonances and in Table 2 residues at resonance 
poles and in Table 3 the coefficients cf^^'^ {i ~ 1,2; j = 2,3,4; / = 0, 1) in the expansion in 
terms of the effective coupling constant an of QCD defined by ([55]) . The number of flavor is taken 
as n/ = 3. Imi^^ is cut-off at A§ = 0.779 GeV^ and the Breit-Wigner formulas at Ai = 26.0 GeV. 
The QCD parameter is fixed at A = 0.216 GeV. The other parameters are determined as follows: 
Case I: to = 4//^ h = 0.243x10^ GeV^, t2 = 0.237x10^ GeV^, Ues = 0.2260x10^ GeV^, tg = 
0.202 xlO^ GeV^. 

Case II: The same as in the case I except for tres = 0.2253x10^ GeV^. 

The value of is obtained to be Xtot = 393.4 for the case I and Xtot — 308.7 for the case II, 
which includes both the data of space-like and time-like regions. The total number of data is 245 
for the Case I and 236 for the Case II. Number of parameters is 36 so that DOF/xmm 1-88 for 
the Case I and 1.54 for the Case II. 
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Table 1: Masses and widths determined by the analysis for the cases I and II. 





case I 


case II 






mass 


width 


mass 


width 


isospin 


n 


(GcV/c2) 


(GcV) 


(GcV/c2) 


(GcV) 




1 


1.341 


0.3221 


1.352 


0.325 




2 


1.379 


0.2204 


1.370 


0.220 


/= 1 


3 


1.599 


0.2636 


1.587 


0.264 




4 


1.824 


0.3679 


1.826 


0.368 




5 


2.048 


0.3848 


2.100 


0.398 




1 


0.78256 


0.844x10-^ 


0.78256 


0.844x10-2 




2 


1.01945 


0.426x10^2 


1.01945 


0.426x10-2 


/ = 


3 


1.212 


0.1582 


1.206 


0.1584 




4 


1.437 


0.2102 


1.440 


0.2104 




5 


1.505 


0.1281 


1.510 


0.1285 



Table 2: The coefficients a['", residues at the resonance poles, determined by the analysis for 



the cases I and II. 





case I 


case II 


isospin 


n 


a('"(GeV2) 


a^'"(GeV2) 


a{'"(GcV2) 


«2'"(GcV2) 




1 


-4.66 


8.45 


-4.47 


8.37 




2 


8.489277 


-17.50252 


7.4739945 


-15.8900 


1=1 


3 


-9.623356 


13.46278 


-8.050218 


10.93951 




4 


7.065036 


-7.310033 


6.091668 


-6.071918 




5 


-0.140 


1.36 


-0.118 


1.11 




1 


0.899887 


0.02568286 


0.8762127 


0.09468219 




2 


-3.625433 


0.5913331 


-3.514823 


0.3151743 


1 = 


3 


7.385961 


-2.033127 


6.954618 


-1.526721 




4 


-3.934473 


-1.019970 


-3.579443 


-2.125879 




5 


-1.028184 


2.582630 


-1.038278 


3.394527 



Tabic 3: The coefiicients cff^'^ of the QCD terms for the cases I and II determined by the x^ 

analysis. 





(:as(> I 


isospin 


i 


ocn.i 

i.2 


QCD,l 


QCDJ 


/ = 1 


1 


0.5505731 


-4.12 


-6.50 




2 


3.758361 


-0.4002 xl02 


0.6224x102 


/ = 


1 


1.108707 


-2.76 


-0.7045 xl02 




2 


-5.215940 


0.6706x102 


-0.19908x10^ 




case II 


isospin 


i 


QCU,I 


QCD,! 


QCU,1 


1=1 


1 


-0.7186148x10-1 


1.48 


-6.99 




2 


4.252983 


-0.4375 xl02 


0.5543 xl02 


1 = 


1 


0.8918455 


-1.10 


-0.6787000 xl02 




2 


-5.617625 


0.7029x102 


-0.19551x10^ 



9 




gl I 

10"^ 10"^ 10° 10^ 

[(GeV/c)^] 

Figure 1: Proton magnetic form factor for the space-like momentum. The sohd curve is the result 
for case I and the dashed one for the case II. 

We illustrate in Figs. 1-9 the calculated results for the form factors. The results for the 
Case I is given by the solid curve and Case II by the dashed one. Figs. 1-4 the results for the 
space-like momentum are illustrated: Fig. 1 the proton magnetic form factors G^j/fipGo, Fig. 2 
proton electric form factor G^/Gd, Fig. 3 the neutron magnetic form factor G^///i„ and Fig. 4 the 
neutron electric form factor. In Fig. 5 we illustrate the ratio of proton electric and proton magnetic 
form factors ^pG^/G^^. We find that G| = at = 6.57 (GcV/c)^ for the case I and = 6.79 
(GeV/c)^ for the case II. The form factor for the time-like momentum |G| is given in Fig. 6 for 
the proton and in Fig. 7 for the neutron. The result for the proton form factor agrees with the 
experimental data, but for the neutron the calculated one becomes larger than the experiments for 
large Q^. 

In Fig. 8 we compare the calculated result for the neutron magnetic form factor G\J ^uGd 
with the recent experiments. The solid curve agrees with the experimental data very well. The 
dashed one becomes a little larger than the result obtained by the CLAS collaboration. However, 
the deviation is not very large. In Fig. 9 we illustrate the result for \G^^/^XpG\j\ for the time- 
like momentum. There seems to be some discrepancy between the experimental data: The ratio 
obtained by Bardin el al. [8] is smaller than that of Aubert et al. [9]. Our result coincides with 
the result of Bardin et al. for small and that of Aubert el al. for large Q^. 

5 Concluding Remarks 

The experimental data for the neutron magnetic form factor for the space-like momentum with 
= 1.4 — 4.8 (GeV/c)^ [7], mentioned in Sec. 1, are reproduced very well by our calculation. 
The absorptive parts of the form factors for the asymptotic region are approximated by the 
power series in the effective coupling constant of QCD, which begins 0(a|j) as is given in 
We have taken three terms in the expansions; the terms of order up to 0{a%) are necessary to 
reproduce the experiments as in the case of deep inelastic electron scattering processes. 
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Figure 2: Proton electric form factor for the space-like momentum. The solid curve is the result 
for case I and the dashed one for the case II. 
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Figure 3: Neutron magnetic form factors for the space-like momentum. The solid curve is the 
result for case I and the dashed one for the case II. 
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Figure 4: Neutron electric form factor for the space-like momentum. The solid curve is the result 
for case I and the dashed one for the case II. 




Figure 5: Ratio of the electric and magnetic form factors of proton for the space-like momentum. 
The solid curve is the result for case I and the dashed one for the case II. 
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Figure 6: Proton form factors for the time-like momentum. The soUd curve is the result for case I 
and the dashed one for the case II. 
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Figure 7: The neutron form factor for the time-like momentum. The solid curve is the result for 

case I and the dashed one for the case II. 
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Figure 8: Neutron magnetic form factor for the space-like momentum in the few (GeV/c)^ region. 
The soUd curve is the result for case I and the dashed one for the case II. 




Figure 9: |G^//UpG^| for the time-like momentum. The solid curve is the result for case I and the 
dashed one for the case II. 
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It is remarked here that the electromagnetic form factor of bosons, both for the space- hke and 
timc-hke momentums, can be explained with recourse to the superconvergent dispersion relation 
with the QCD constraints [TS] . 

For the electric form factor of proton there are deviation of the dispersion theoretical calculation 
from the experimental data for large Q^, where the data were obtained by using the Rosenbluth 
formula. The discrepancy may imply the necessity of correction of two photon processes to the 
experimental data |3T] |42) . 

We used the experimental data for the helicity amplitudes obtained by Hoher and Schopper in 
which the contribution from the p meson is included. As their data are limited to low t (< 0.779 
(GeV/c)^), we do not have sufficient data for the region s < 4m^. We supplemented the unphysical 
region for 1=1 state by introducing vector bosons with the mass, my < 1.4 GeV/c^. For the 
isoscalar state we also introduced a vector boson with the mass about 1.2 GeV/c^. 

In our calculation we treated all of the vector boson masses and widths as parameters. If 
they are kept at experimental values, we get poor results. The superconvergence conditions are so 
strong that the value of is very sensitive to the mass and width. The masses are obtained to 
be smaller than the experimental value and the existence of vector bosons with the masses around 
1.2 ~ 1.4 GeV/c^ are necessary both for the 1=1 and / = states. 

To conclude the paper we remark on the mass around 1.2 GeV/c^. We have introduced the 
vector boson to supplement the lack of information on the the small Q^. However, both for / = 
and 1=1 states there are indications of resonances observed by the processes e'^e" — rjir'^n^ , 
7P —7- LWK^p and B — D*un:~ [43j . Incorporation of further resonances may improve results for 
the time-like momentum. 
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